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Abstract 

In this paper we propose the use of neural interference as the origin of quantum- 
like effects in the brain. We do so by using a neural oscillator model consistent 
with neurophysiological data. The model used was shown to reproduce well 
the predictions of stimulus-response theory. The quantum-like effects are pro- 
duced by the spreading activation of incompatible oscillators, leading to an 
interference-like effect mediated by inhibitory and excitatory synapses. 
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1. Introduction 



Quantum mechanics is one of the most successful scientific theories in his- 
tory. From it, detailed predictions and extremely precise descriptions of a wide 
range of physical systems are given. Thus, it is not surprising that researchers 
from other disciplines started asking whether the apparatus of quantum mechan- 
ics could be used t o help with the description of social or behavioral systems 
( Khrennikovll2010l) . [n applying the quantum formalism, most researchers do 



not claim to have found social or behaviora l phenomena that are determined 



by the physical laws of quantum mechanics (jBruza et al. . 120091) . Instead, they 



claim that by using the representation of states in a Hilbert space, together 
with the corresponding non-Kolmogorovian interference of probabilities, a bet- 
ter description of observed phenomena could be achieved. To distinguish the 
idea that the formalism describing a system is quantum without requiring the 
system itself to be, researchers often use the term quantum-like. 

One of t he first applicat i ons of quantu m-like dynamics out side of physics was 
in finances (iBaaq uiel 1997 ; Haven . 20021) and in p sychology (jBusemever et al "" 



l2006t lAertsl l2009f h In their work, iBaaauid (|l997f ) showed that a more general 
form of the Black- Scholes option pricing equation was given by a Schrodinger- 
type equation, and iHaven (2002| linked the equivalent of the Plank constant to 
the existence of arbitrage in the model. In psychology, most of the initial ef- 
fort was focused on using quantum-like dynamics and Hilbert-space formalisms 
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to describe discrepancies between data and the " c lassica l " models using stan - 



dard probability theor y (see, for example, Aerts 120091); Asano et al 



Busemcv er et al.l (12006T) . iBusemcvcr and Wand (12 007): Busem ever et al 



2010); 



2009); 



Khrennikovl (|2007l . l2009l) ; iKhrennikov and Havenl (|200^ ; iKhrennikovl (|201lK 



To try to descri be quantum-like effects in cognitive sciences, in a recent paper 



Khrennikovl (|201lh used a classical electromagnetic-field model for the mental 
processing of information. In it, Khrennikov argues that the collective activity 
of neurons involved in mental processing produces electromagnetic signals. Such 
signals the propagate in the brain, leading to interference and the appearance of 
quantum-like correlations. He then uses his model to show how the quantum- like 
effects could help in understanding the binding problem. 

In this paper we present an alternative model exhibiting many of the char- 
acteristics of Khrennikov's, such as signaling, superposition, and interference 
but that does not rely on the propagation of electromagnetic fields in the brain. 
Instead, we show that those effects could be obtained by directly considering the 
activities of collections of firing neurons coupled via inhibitory and excitatory 
synapses. The model we present was initially proposed as an attempt to obtain 
behavioral stimulus- response theory from the neuronal activities in the brain 
( Suppes et al.l . 120121) . Here, we show that by using it we obtain quantum- like 



effects, and we use it to model a violation of Savage's sure-thing principle. 

This paper is o rganiz ed as follows. In section [5] we introduce Suppes, de 
Barros, and Oas 's (|2012h oscillator model, with special emphasis to its neuro- 
physiological interpretation and connection to SR theory. Then, in section [3J 
we present a contextual t heory of quantum-like pheno mena coming from classi- 
cal interference of waves ( de Barros and Suppes . [2009h . Finally, in section [2] we 
show computer simulations of neural oscillators and how they fit some quantum- 
like experimental data in the literat ure. We end with some remarks about the 
relationship between our model and IKhrennikov ((2orl . 



2. SR-theory neural oscillators 



In this section, we start by briefly reviewing S R theory, and then we present 
the oscillator model used bv lSuppes et al.l (|2012r ) to model it. Our presentation 
is not intended to be complete; for example, we do not discuss the stability 
of solutions or the dynamics learning. Instead, we focus on the main features 
of the model that are relevant for the e mergence of quantum -like effects. The 



interested reader is referred to reference (jSuppes et al.l . l2012f ) for details 



In m a them atical psychology, one of the most successful theories of learning 
is Estes (1950) stimulus-response (SR) theory, and we chose to model it for 
several reasons. First, it fits extremely well a multitude of behavioral exper- 
iments. Second, it makes definite and, to some, surprising predictions, such 
as probability matching and asymmet r ic con ditional probabilities for symmet- 



ric reinforcements (see ISuppes et al.l . l2012l and references therein). Third, 
SR theory is built on a simple an d well-defined trial structure, which allows 
for a formal axiomatic treatment ( Suppesl . 20021 ). This formal mathematical 
structure provides the necessary tools for the proof of important representation 



2 



theorems. Finally, despite its simplicity and rigid trial structure, it is possible 
to prove that SR theory is powerful enough to pr ovide represen tations of finite 
automata and, consequently, of natural language riSuppesl . l2002h . Therefore, SR 
theory offers a framework that is both formally simple and yet powerful enough 
to account for many higher-level cognitive processes. 

SR theory, in its continuum of respons es version, can be fo rmalized in terms 



of a stochastic process (here we follow ISuppes et all 120121) . Let (fi, J-, P) 
be a probability space, and let Z, S, X, and Y be random variables, with 
Z : fi -> £l 5 l S : n -> S, X : SI -> R, and Y : H -> E, where S is the set of 
stimuli, R the set of responses, and E the set of reinforcements. Then a trial in 
SR theory has the following structure. 



Z n — > S„ — > X n — > Y„ 



(1) 



>; 



Intuitively, the trial structure works the following way. Trial n starts with a cer- 
tain state of conditioning, represented by the random variable Z„ = (z^ 

The vector ^j"',...,zt'J associates to each stimuli Sj € S, i = 1, 
where m — \S\ is the cardinality of S, a value 



. .,771, 



» 



Si is sampled with probability P (S., 



on trial n. Once a stimulus 
Si\sieS) = — , its corresponding 



z- determines the probability of responses in R by the probability distribution 



dx, 



K (r|4 n) ) , i.e. P ( 01 < X„ < a 2 |S„ = s h Z n4 = z^) = Q k 

where k \ x\zf^ is the probability density associate to the distribution, and 

where Z„.i is the i-th component of the vector (jz^, . . . , zfj . The probability 

distribution K \ r\z^j is the smearing distribution, and it is determined by its 

variance and mode z| n . The next step is the reinforcement Y„, which is ef- 
fective with probability 8, i.e. P ^Z„ + i.i = y\S n = S{, Y„ = y, Z n/ 



y (n) 



and P ^Z n +i.j =zj^|S n =Sj,Y„ = y, Z nj j = z| n ^ = 1 — c. The trial ends with 

a new (with probability c) state of conditioning Z ra _)_i . 

Because of its simplicity and focus on behavioral outcomes, mathematical 
SR theory does not provide a clear connection to t he pro cessing of information 
by the brain. To bridge this gap, ISuppes et al. ( 2012 ) proposed a response 
computation model inspired by SR theory. In their model, they provided a 
neural processing model able to reproduce the main stochastic features of SR 
theory, including the conditio nal probabilit i es for a continuum of responses. 

The main characteristic of lSuppes et al] (|2012l ) is the use of neural oscillators 
for computing responses after learning. There are many differ ent ways to use 
neural oscillators to model lear ning (see Vassiliev a et all I201U and references 
therein), but to our knowledge. ISuppes et al.l (|2012l ) is the only one that makes 
a clear connection between neural oscillator computations and the behavioral 
responses. Furthermore, their model is simple enough to help understand the 
basic principles behind the computation of a response from a stimulus without 
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compromising the neurophysiological interpretation of it. This must be empha- 
sized, as most neural models constructed from a neurophysiological basis are so 
computationally intensive and complex that it is very hard to understand why 
or how one response was select i nstead of another. A s we shall see later in this 
section, this is not the case with lSuppes et al.l (|2012i h 

Let us start the description with the motivations for choosing the particular 
setup used. The interest is to model higher cognitive functions, such as lan- 
guage or decision making. In language, there is a large body of work suggesting 
that syllables and words are represented on electroencephalogram (EEG) wave 
patterns (see Supped . 2002, Chapter 8, and corresponding references). EEG 
waves, on the other hand, are the res ult of a large assembly of ne urons firing 
synchronously or quasi-synchronously ( Nunez and Srinivasanl [2006) . Thus, we 
can think of, say, a word as represented in the brain by a set of coupled neurons, 
such that when activated by an external stimulus they fire coherently. Such an 
ensemble of synchronizing neurons could, in first approximation, be represented 
by a periodic function Fr(t), where T is the period of the function^. 

One of the interesting things about an ensemble of synchronizing neurons is 
that, when coupled to another ensemble, even if weakly, they may synchronize. 
To understand how this happens, let us look at the qualitative behavior of 
individual neurons in two ensembles, A and B, each with a large number of 
neurons firing coherently. Let ha be a single neuron in ensemble A firing with 
period T, and let tib be another neuron in B firing with period T". Let t' be a 
particular time such that ns would fire if it were not connected to ha (i-e., the 
time it would fire because of its natural periodicity T"). Let us also assume that 
ha is connected to tig via an excitatory synapse. Because of this excitatory 
connection, the firing of ua at a time t j§ t 1 would have the effect to anticipate 
ns's firing to a time t" closer to t, thus approaching the timing of firing for both 
neurons. If more neurons from ensemble A are coupled to ns, the stronger the 
effect of anticipating its firing time. Furthermore, since ub is also coupled to 
neurons in A, its firing has the same effect on A (though this effect does not need 
to be symmetric). So, the couplings between A and B lead both periods T and 
T' to approach each other, i.e., A and B synchronize. In fact, it is possible to 
prove mathematically that, under the right conditions, if the number of neurons 
is large enough, the sum of the several weak synaptic interactions can c ause 
a strong effect, making all neurons fire close together (|lzhikevichL l2007t ). In 
other words, even when weakly coupled, two neural ensembles represented by 
oscillators may synchronize. 

We now turn to a more rigorous mathematical description. Our underlying 
assumption is that we have two neural ensembles A and B described by periodic 
functions Fx A (t) and Ft b (t) with periods Ta and Tg, respectively. Without 
loss of generality in our argument, we assume that Fj- A (t) and Fj- B (t) have 
the same shape, i.e., there exists a constant c such that Fx A (t) = Fx B (tc). 



1 We remark that neurons in this ensemble do not need to fire at the same time, but only 
that their frequencies synchronize. They may still be in synchrony but out of phase. 
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A simple way to describe the synchronization of A and B is to rewrite the 
time argument in the functions, rewriting them as Ft a (tcxA) an d Ft b (tas)- 
Clearly if no couplings exist between A and B, their dynamics is not modified 
and a a = «b = 1. However, when the ensembles are couples, a a and olb 
are a function of time, their synchronization means that they evolve such that 
Ft a (totA (t)) — Ft b {tctB (£))■ Thus, when we are studying the dynamics of 
coupled neural oscillators, it suffices to study the dynamics of the arguments of 
Ft a and Ft b above, which we call their phases. 

To illustrate this point, let us imagine a simple case where when there is no 
interaction A follows a harmonic oscillator with angular frequency lua — 2it/Ta- 
Then Ft a (t) — A cos [u)At + 5a), where and 5a is a constant, and expressions 
hold for Ft b (t). Following the above paragraph, we rewrite 



where 



F TA (t)=Aeos(<p A (*)), 



ip a (t) = oj A t + 5 A (2) 



is the phase. Let us now focus on the p hase dynamics. 

For a set of coupled phase oscillators, Kuramo to ( 1984) proposed a simple set 



of dynamical equations for ifA (t)- To better understand them, let us assume 
that in the absence of interactions if a evolves according to @. Thus, the 
differential equation describing tpA would be 



dp a (t) 
dt 

and for ipg, 

dtp B (t) 



^A, (3) 



uj b . (4) 



dt 

In the presence of a synchronizing interaction, Kuramoto proposed that equa- 
tions © and (01 should be replaced by 



and for tps, 



d(fiA (t) 
dt 

dip B (t) 



^a - k A B sin [ifA (t) - (Pb (*)) , (5) 



ub ~ k B A sin (tp B (t) - ip A (t)) . (6) 



dt 

It is easy to see how equations and ^ work. When the phases cfA and ipe 
are close to each other, we can linearize the sine term, and rewrite equations 
© and © as 

dtp A (t) 



and for ips, 



dt 

dipB jt) 
dt 



w A t - k A B {<Pa (t) - Pb (*)) , (7) 
io B t - k B A (pb (t) - pa (*)) • (8) 



The effect of the added term in equation is to make the instantaneous fre- 
quency increase above its natural frequency ua if <Pb is greater than if a-, and 
decrease otherwise. In other words, it pushes the phases if a an d ifB towards 
synchronization. We should emphasize that Kuramoto's equations were pro- 
posed for a large number of oscillators, and his goal was to find exact solutions 
for such large systems. However, as we described above, his equations also 
make sense for systems with small number of oscillators (for applications of 
K uramoto's equ a tions to systems with small numb e r of o s cillato r s, in addition 



to ISuppes et all (120121) see also iBillock and Tsoul (|2005l 1201 if ); ISeliger et af. 
( 20021) : iTrevisan et al.l (|2005t ) and references therein) . For a set of N coupled 



phase oscillators, Kuramoto's equations are 



= u>j - kij sin (fi (t) - tfj (*)) . (9) 



dt 

3=1 

3+i 



But independent of whether we have large or small number of oscillators or not, 
it is possible to prove that for oscillating dynamical systems near a bifurcation 
(which is the case for many neuronal models), Kuramoto's equatio ns present a 
first order approximation for the synchronizing coupled dynamics ( Izhikevichl 
l2007h . 



In the model proposed in (jSuppes et al. . 2012 ). each stimulus in the set S of 



stimuli corresponds to an oscillator, i.e. there are m stimulus neural oscillators, 
{si (£)}, i — 1, . . . , N. Here we use the notation Sj (t) to distinguish between the 
oscillator and the actual stimulus s, G S. Once a distal stimulus is presented, 
the processing of the proximal stimulus leads, after sensory processing, to the 
activation of neurons in the brain corresponding to a neural oscillator repre- 
sentation of such stimulus. Once a stimulus oscillator is activated by starting 
to fire synchronously, a set of response oscillators, r\ (t) and (t), connected 
to the active stimulus oscillator via the couplings k Si _ rj and fc ri ,r 2 lead to the 
synchronization of the stimulus and response oscillators. The couplings k Si , rj 
and £Vi,r 2 correspond to the state of conditioning in SR-theory. Depending on 
the details of the synchronization dynamics, as explained below, a response is 
selected (this being the equivalent of sampling X n ). Finally, during reinforce- 
ment, a reinforcement oscillator e y (t) is activated, and its couplings with Si (t), 
ri (£), and ri (£), together with k ai)Tj and k ri) r 2 , leads to a dynamics that allow 
for changes in k Si >r . and k ri <r2 , which accounts for the last step in the SR trial 
([1} , with a new state of conditioning. 

Before we detail the dynamics of the oscillator model, it is important to 
discuss how Sj (i), r% (t), and T2 (t) can encode a response through their cou- 
plings, and how can we interpret such couplings. Let us write the case where 
Si is sampled, and let us assume that once the dynamics is acting, all oscilla- 
tors synchronize, acquiring the same frequency (though perhaps with a phase 
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difference). Then, 

Si (t) = Acos (tp Si (i)) = Acos(uj a t), (10) 
ri (t) = Acos((p ri (t)) = Acos(wo* + ^i), (11) 
r 2 (t) = Acos(yv 2 (f)) = Acos(w t + ^2), (12) 

where s(t), ri(i), and rzit) represent harmonic oscillations, ip s (t), tp ri (t), and 
ifr 2 (t) their phases, 5<f>i and are constants, an d A their amplitude, assumed 
to be the same. As ar gued in (ISuppes et all 20121). since neural oscillators have 



a wave-like behavior (|Nunez and Srinivasanl 2006), their dynamics satisfy the 



(f(t)) t = Jim - / f(t)dt 



principle of superposition, thus making oscillators prone to interference effects. 
As such, the mean intensity in time gives us a measure of the excitation of the 
neurons in the neural oscillators. For the superposition of s(t) and Ti(t), 

h = ((*(*) +»i(*)) a ) t 

= ( S (t) 2 ) t + (r 1 (tr) t + (2s(t)r 1 (t)) t , 

where 

1 ' T 

T^o T 

is the time average. It is easy to compute that 

h = A 2 (1 + cos (5<j>i)) , 

and 

h = A 2 (l + cos(^ 2 )). 

From the above equations, the maximum intensity for each superposition is 2 A 2 , 
and the minimum is zero. Thus, the maximum difference between I± and I2 
happens when their relative phases are tt. We also expect a maximum contrast 
between I\ and I2 when a response is not in between the responses represented 
by the oscillators r\ (t) and r 2 it). For in between responses, we should expect 
less contrast, with the minimum contrast happening when the response lies on 
the mid-point of the continuum between the responses associated to r\(t) and 
r 2 (i). The ideal balance of responses happen if we impose 

Sfa = 5<f>2 + tt EE 5<f>, (13) 

which results in 

I x = A 2 (1 + cos ((50)), (14) 

and 

I 2 = A 2 (1 - cos ((50)) . (15) 

From equations (fT4|) and (fl~5|) , let x € [— 1, 1] be the normalized difference in 
intensities between n and r 2 , i.e. 

X EE y 12 = COS (Slfi) , (16) 
il + 1 2 
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< Sip < 7r. The parameter x codifies the dispute between response oscillators 
r\ (t) and it). So, we can use arbitrary phase differences between oscillators 
to code for a continuum of responses between —1 and 1. 

As we mentioned above, the dynamics of the phase oscillators leading to 
(p~0t> (fT2]) are modeled here by Kuramoto equations ©. However, as we also 
mentioned, (0) leads to the synchronization of all oscillators with the same 
phase, which is not what we need to code responses according to equation (|16p . 
Naturally, Kuramoto's equations need to be modified to encode the phase dif- 
ferences Sip in pU)) - (fT^l) . This can be accomplished by adding a term inside the 
sine, resulting in 

dip-it) N 

— — = cjj - hj sin fa (t) - ipj (t) + Sij) . (17) 

The main problem with such change is how to interpret it. Equation © had a 
clear interpretation: the firing of oscillator neurons coupled through excitatory 
synapses displaced the phase of the other oscillator into synchrony. But this 
interpretation does not make sense for (fTT)) . Why would oscillators be brought 
close to each other, but be kept at a phase distance of Sij? 
To understand the origin of Sij , let us rewrite (fT7|) as 

JV 



dt 



} J hj cos (Sij) sin (tpi - ipj) 



A' 



- hj sin (S^) cos ((fi - ipj) . (18) 

Since the terms involving the phase differences Sij are constant, we can write 
p| as 

= <^ - ^ [kfj sin (^ - + kjj cos (^ - ipj)) , (19) 
i=i 

where kfj = hj cos (Sij) and /c/^ = A;^ sin (Sij). Equation (|17[) now has a clear 
interpretation: fcy makes oscillators ip t (t) and ipj (t) approach each other, kf- 
makes them move further apart. Thus, kfj corresponds to excitatory couplings 
between neurons and kfj to inhibitory couplings. In other words, we give mean- 
ing to equations (| 1 7[) by rethinking their couplings in terms of excitatory and 
inhibitory neuronal connections. 

To summarize the response process, once a stimulus oscillator s,; (t) is acti- 
vated at time t an on trial n, the response oscillators are also activated. Because 
of the excitatory and inhibitory couplings between stimulus and response os- 
cillators, after a certain time their dynamics lead to their synchronization, but 
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with phase differences given by (fTTJ)) . Here we point out that, due to stochas- 
tic variations of biological origin, the initial conditions Sj (t StU ), r\ (£ s ,n)j an d 
f2 (t s ,n) vary according to the distribution 

Therefore, the phase differences at the time of response, t Xjn , may not be exactly 
the ones given in (|13[) . which gives an underlying explanation for a component 
of the smearing distribution^. Interference between oscillators, determined by 
their phase differences, leads to a relative intensity that codes a continuum 
of responses according to x in equation (|16p . The variable x is the response 
computed by the oscillators. 

We now turn to learning. Here we will only present the main features of 
learnin g relevant to this pa per. The actual details are numerous, and we refer 



the to (jSuppes et al.l . |2012| ). For learning to happen, the couplings kf^ and k^ 



need to change during reinforcement. Since reinforcement must have a brain 
representation in terms of synchronously firing neurons, during reinforcement 
a neural oscillator <p y (t) gets activated. We assume that (p y (t) is a fixed rep- 
resentation in the brain, with frequencies that are independent of the stimulus 
and response oscillators, and therefore can be represented by 

(fy (t) = UJ e t. 

Furthermore, since ip y is coding a response y being reinforced, from (|16p the 
phase of the oscillators must be related to y via 

Sip = arccosy. 

During the reinforcement period, the (active) stimulus and response oscillators 
evolve according to the same equations as before, but now with an extra term 
due to the reinforcement oscillator, according to 

-J^ = uii - ^2 k ij sin ~ Vi) 
- ^2 k E sin (<pi - tpj) 

Si,ri,r 2 

-K Q sin (tp Sj - uj E t + Ai) , (21) 



2 We note that we here only model the brain computation of a stimulus and response, and 
not the representation of the distal stimulus into brain oscillators and the representation of 
responses into actual muscle movement or sound production. Of course, those extra steps 
between the stimulus and the response will lead to further smearing of the actual observed 
distribution, and a detailed brain computation o f K f r|z| n ^ J should include such factors. The 

interested reader should refer to ISuppes et al.l l|2012T l for details about the used parameters 
and the fitting of the model to experimental data. 
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where A y - reflects the reinforced response y, 



Aj = - (<5i, ri + <5 ljr2 ) arccosy + ^, r2 7r, 
and where is Kronecker's delta. During reinforcement, c ouplings also chang e 



according to a Hebbian rule represented by the equations (Se liger et al.l . 1200 



dk E - 

-tt = e ( K o) [a cos fa -<fj)- hi 



(22) 



and 



— = e{K ) [a sin ((fi-^j) 



(23) 



where 



0iiK o < K' 
6q otherwise, 



(24) 



where eo <C ujq , a and Kq are constant during At K , and K' is a threshold 



constant throughout trials ( Hoppensteadt and Izhikevich , 1996 b. a). We also 
assume that that there is a normal probability distribution governing the cou- 
pling strength Ko between the reinforcement and the other active oscillators. It 
has mean Kq and standard deviation <jk - Its density function is: 



f(K ) = 



'2-K 



exp 



2< 



(25) 



So, when a reinforcement oscillator is activated, the system evolves under the 
coupled set of differential equations (|2"TT) . (f2"2"j) . and ([2"B")) . 

To summarize, we described a model of SR theory in terms of neural oscilla- 
tors. In such a model, stimulus and responses are represented by collections of 
neurons firing in synchrony. Such collections of neurons, when coupled, synchro- 
nize to each other. The relative phases of synchronization lead to interference 
of sets of neurons, and such interference results in the coding of a continuum of 
responses. We left our many of the details of the stochastic processes involved 
in this mode l, as they are not ne cessary for this paper, but details can be found 
in reference (jSuppes et all 120121 ). 



3. What are quantum-like effects? 

In this section we discuss possible quantum-like effects in the brain and 
how they can be understood in terms of oscillators. Before we proceed, it is 
important to discuss what we mean by quantum-like effects. 

Quantum mechanics is one of the most successful theories in history, ac- 
counting for numerous phenomena in nature. However, such success did not 
come without a price, as quantum mechanics presents a view of the world that 
most people would find disturbing, to say the least. This was certainly the view 
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of many of the founders of quantum mechanics, such as Plank, Einstein, and 
Pauli. At the core of their concerns were three characteristics of quantum me- 
chanics: nondeterminism, contextuallity, and nonlocality. Let us examine each 
one of those characteristics, and see how they can be relevant to brain processes. 

We start with nondeterminism. Intuitively, a dynamical system is determin- 
istic if the current state of the system completely determines the future state of 
such system, and it is nondeterministic otherwise. Of course, for practical pur- 
poses, there are many nondeterministic processes, but physicists used to believe 
that the underlying dynamics for such processes could be deterministic. For ex- 
ample, when we toss a coin or throw a die, physicists believe that the newtonian 
dynamics would allow for the computation of the final state of the die or coin, 
if we were to know exactly the initial conditions, i.e., the state of the system. 
Whether we can distingu ish between a determin i stic a nd stochastic dynamics is 
not the important issue (jSuppes and de Barrosl I1996I) . The main point is that 
for quantum theory the underlying dynamics is essentially nondeterministic. 

It should be clear that when we talk about brain processes, we are not 
concerned about quantum-like nondeterminis tic effects. Such effects are ofte n 
modeled in psychology stochastic processes ( Busemever and Diederich . 2010h . 



For example, the stochastic SR theory presented above has a nondeterministic 
characteristic, yet it is not concerned with the underlying sources of nonde- 
terminism. Furthermore, suc h stochasti c proc esses may come from inhering 



biological noises in the brain ( Josic et al. . 20091 ). Insofar as nondeterminism is 



concerned, there is no need to add further quantum-mechanical mathematical 
or conceptual machinery to account for it in the brain. 

Now let us examine nonlocallity. To discuss nonlocallity, we need to discuss 
causality. Let A and B be two events, and let A and B be corresponding ±1- 
random variables corresponding to whether A and B occur red. We say th at C is 
a prima facie cause of A if P(A = 1|C = 1) > P (A = 1) ( Suppes . 1970l) . Now, 
it is possible to prove that for some quantum systems, there are space- like events 
A and A' such that A is a prima facie cause of A' . What is key in the previous 
statement is that A and A' are space-like separated events; in other words, 
whatever mechanism is making A' influence A must be superlumina{f|. This, it 
may be argued, is the main issue with quantum processes: to understand them 
we need to accept superluminal causation. However, we believe that quantum- 
like effects in the brain are not the outcome of superluminal causation. The 
brain is small, of the order of tens of centimeters, and an electromagnetic field 
would only need of the order of 10~ 10 seconds to travel the whole brain. Given 
that brain processes are orders of magnitude slower than 10~ 10 seconds, no 
phenomena in the brain could not be accounted for non-superluminal causality. 
In other words, we doubt it would be possible to create an experiment where 
nonlocal effects in the brain would be detected in a superluminal way. 



3 We are summarizing a very complex and subtle discussion in a few sentences. Th e reader 
intere sted in more details and in a mathematically rigorous treatment is directed to ISuppesI 
i2002ft and references therein. 



11 



Here we need to differentiate our use of nonl ocality from what is some- 



times found in the literature on quantum cognition ( Nelson et al. . 20031). Some- 



times , authors consider a violation of some form of Bell's inequalities (Belli. 11964 , 
19661) as evidence of nonlocality. However, from our discussion, the correlated 
quantities must be measured such that no superluminal signal explaining the 
correlations is possi ble, otherwise a viola tion of Bell's inequalities would only 



imply contextuality (jSuppes et all Il9 96a). Classical fields, such as the electro- 



magnetic fie ld used by Khrenniko v, also exhibit contextuality, violating Bell's 
inequalities (jSuppes et al. . Il996bl) . Here we are using nonlocality in the strict 
sense of "spooky" correlations that cannot be explained without using superlu- 
minal interactions. 

The last quantum-like issue we discuss is contextuality. Early on, the wave 
description of a particle led physicists to realize, through Fourier's theorem, that 
it was impossible to simultaneously assign values of momentum and position 
to it. This developed into the concept of complementarity. Two observables 
properties are complementary when they cannot be simultaneously observed, 
i.e., when the corresponding hermitian operators 0\ and O2 for such properties 
does not commute ([Oi,C>2] 7^ 0). But the question is not whether we can 
measure 0\ and O2 simultaneously, but whether we can assign those properties 
values even when we cannot measure them. In other words, can we make a 
table that fills out the values of 0\ and O2 for every trial? If the number 
of observab les is large enou gh, the answer to that que stion was given in the 
negative by iBell ( 19661 ) and Kochen and Specker ( 1975 ) , who showed that any 
attempts to fill out such tables would result in values inconsistent with those 
predicted by quantum mechanics. In other words, there are quantum systems 
whose observables cannot be assigned values when we do not measure them, 
and therefore we cannot provide a consi s tent joint p r obabi l ity distribution fo r 



and tncretore we cann ot provide a consi s tent lomt p r obabi l ity distribution t or 
such observables (jde Barros and Supped . I2OO0I |200ll l20ld: Ide Barrosl l201ll Fl 



This impossibility of defining a joint probability distribution is related to the 
change of values of the random variables when contexts change. For example, a 
random variable Oi representing 0\ the outcomes of an experiment may not be 
the same random variable when O2 is being measured. Such random variables 
change with the context. Contextuality is not new in physics, nor in psychology. 
It happens in classical s ystem s, such as with the classical electromagnetic field 
(|de Barros and Suppe 1 120091 ). However, we claim it is contextuallity, with its 
associated impossibility of defining joint probability distributions, that leads to 
quantum-like events in brain processing. 

At this point, it is important to raise the issue of how much quantum me- 
chanics apparatus we may need to describe brain processes. Quantum mechanics 
brings to the table more than just contextuality, nonlocality, and nondetermin- 
ism. It also brings an additional mathematical structure, formalized by the 
algebra of observables on a Hilbert space. To see this, let us examine the fol- 



4 Thou gh with the redefinitio n of what constitutes a single photon, local models are possible 
(see, e.g. ISuppes et all [l996cl) . 
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lowing case, analyzed in details by ISuppes and Zanotti ( 1981 ); Suppes et al 



(1996a). Let X, Y, and Z be three ±l-valued random variables with observed 



pairwise joint expectations given by 

E (XY) = E (XZ) = E (YZ) 



where e > 1/3. ISuppes and Zanottil ( 198ll ) showed that such expectations are 



not consistent with the existence of a joint probability distribution for X, Y, 
and Z if 

-1 < E (XY) + E (XZ) + E (YZ) < 1 + 2 min {E (XY) , E (XZ) , E (YZ)} . 

Thus, if they represent the outcomes of local measurements, they are contextual. 
If, on the other hand, they each represent space-like separated measurement 
events, they are non-local. And each random variable is indistinguishable from a 
tossed coin, and are therefore nondeterministielf|. However, X, Y, and Z cannot 
represent quantum mechanical observables, as, since they all commute, we could 
in principle simultaneously observe X, Y, and Z, which is a sufficient condition 
for the existence of a joint probability distribution^]. But we could, on the 
other hand, imagine a social-science or psychology example (albeit probably a 
contrived one) where corre lations such as th e ones exhibited above could appear. 
Furthermore, as shown bv lde Barrosl (|2012h . such random variable correlations 



can be obtained from the same neural oscillator model from section [2] 

In this section we discussed the different ways in which quantum-like pro- 
cesses may show up in the brain. We argued that among them, the most prob- 
able is contextuality. Contextuality may appear in a system due to several dif- 
ferent reasons. For example, certain systems involve dynamical behavior that 
is dependent on boundary conditions, and those may change with the context. 
Another possibility is when complex systems have different parts that can in- 
terfere, changing the overall behavior of the part s due to changes in other parts. 



The latter possibility is what iKhrennikovl (|2011r ) explored in his field model. In 



the next section, we will show that the complex interaction of neural oscillators 
are also another possible explanation for quantum-like behavior in the brain. 



4. Quantum-like effects for neural oscillators 

The experimental violation of Savage's sure-thing principle (STP) is con- 
sidered an example of a quantum-like decision making, so, let us start with a 
description of the STP. Savage presents the following example. 



5 The distinction b e tween nondeterministic and deterministic is a com plicated one, c.f. 
Suppes and de Barros (1994), and we refrain from it in this paper. See also jWerndll J2009T) . 

6 This comes from the fact that, since they commute with each other, there exists a set of 
orthogonal base vectors in the Hilbert space representation such that the observables corre- 
sponding to X, Y, and Z, namely X, Y, and Z, are all diagonal in this base. 
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"A businessman contemplates buying a certain piece of property. He 
considers the outcome of the next presidential election relevant to the 
attractiveness of the purchase. So, to clarify the matter for himself, 
he asks whether he should buy if he knew that the Republican can- 
didate were going to win, and decides that he would do so. Similarly, 
he considers whether he would buy if he new that the Democratic 
candidate were going to win, and again finds that he would do so. 
Seeing that he would buy in either event, he decides that he should 
buy, even though he does not know which event obtains, or will ob- 
tain, as we would ordinarily say. It is all too seldom that a decision 
can be arrived at on the basis of the principle used by this business- 
man, but, except possibly for the assumption of simple ordering, I 
know of no other extralog ical principle g overning decisions that finds 



such ready acceptance." (jSavagd . Il972l pg. 21) 



The idea illustrated is that if X is preferred to Y under condition A and also 
under condition -^A, then X is preferred to Y, which Savage calls the sure-thing 
principle. 

Putting it in a more formal way, consider the following three propositions, A, 
X, and Y, and let P(X\A), the conditional probability of X given A, represent 
a measure of a rational belief on whether X is true given that the condition A is 
true. Since the measure P requires ra tionality of be l ief, it follows tha t it satisfies 
Kolmogorov's axioms of probability ( Javnesl 2003 : Galavotti . 2005t ). We start 



with the assumption that 

P(X\A) > P(Y\A). 

This expression can be interpreted as stating that if A is true, then X is preferred 
over Y. If we also assume 

P(X\-nA) > P(Y\^A), 
and we obtain, multiplying each inequality by P (A) and P (~>A) , respectively, 

P{X\A)P {A) + P(X\->A)P (->A) > P{Y\A)P (A) + P{Y\^A)P (->A) , 
and from P(Ah^A) = 1 and the definition of conditional probabilities we have 

P{X) > P(Y). 

In other words, if we prefer X from Y when A is true, and we also prefer X from 
Y when —>A is true, then we should prefer X over Y regardless of whether A is 
true. This, of course, is equivalent to the above example when A is equivalent 
to a Democrat president winning, ^A to Republican, and X — —Y is equivalent 
to buying. 



Thoug h STP should hold if agents are making rational decisions. iTverskv and Shafir 



(1992) and I Shafir and Tverskvl fl992) showed that it was violated under certain 



conditions. In their work, they presented several simple decision-making prob- 
lems to Stanford students. For example, in one question students were told 
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about a game of chance, to be played in two steps. In the first step, not volun- 
tary, players had a 50% probability of winning $200 and 50% of loosing $100. 
In the second step, a choice needed to be made: whether to gamble once again 
or not. If a player accepted a second gamble, the same odds and payoffs would 
be at stake. When told that that they won the first bet, 69% of subjects said 
they would gamble again on the second step. When told they lost, 59% of the 
subjects also said they would gamble again. Following the discussion above, we 
translate the experiment into the following propositions. 

A = "Won first bet," 

= "Lost first bet," 

X = "Accept second gamble," 

Y = "Reject second gamble," 



and we have that 
and 



P{X\A) = 0.69 > P(Y\A) = 0.31, 



P(X\^A) = 0.59 > P(Y\-iA) = 0.41, 

since P{X) = 1 — P(Y), as X = -Y. Clearly, X is preferred over Y regardless 
of A. Later on the semester, the same problem was presented, but this time 
students were not told whether the bet was won or not in the first round. In 
other words, they had to make a decision between X and Y without knowing 
whether A or -^A. This time, 64% of students chose to reject the second gamble, 
and 36% chose to accept it. Since in this case 

P{X) = 0.36 < P(Y) = 0.64, 

there is a clear violation of the STP. 

We now turn to t he representation of the ab ove situation in terms of neural 



oscillators. Following lTverskv and Shafirl (|1992n . we call 'Won/Lost Version" the 



case when students were told about A is called and the case when students were 
told that A was hidden from them the "Disjunctive Version." Let us start with 
the Won/Lost Version. In the simplest case, we have two stimulus oscillators 
corresponding to the stimuli associated to the brain representation of A and -<A. 
Let us call s a (t) the oscillator corresponding to A, and s-a (t) to ->A. As before, 
response oscillators are r\ (t) and ri (t). For simplicity, we assume that all 
oscillators have the same natural frequency u>i = ujq. We emphasize that this is 
not a necessary assumption, but since the dynamics will lead to synchronization, 
this will make the overall computations simpler. Then, when s a is activated, so 
are the response oscillators, and the dynamics is given by 



ip Sa = oj - kf a>ri sin(<£ SQ , - (p ri ) 

- k L, ri COs(<yS Sa - tp ri ) 

~ k i a .r 2 COs((^ Sa - tp r2 ) , (26) 
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-kf ltT2 sin((^ ri -<p r2 ) 

-^ iir2 cos(^ n - 99 r3 ), (27) 



¥>r 2 = w - fc^ !ri sin (<y9 r2 - yvj 
-^,»„ sin(^ r2 - ip Sa ) 
-k' r2 ri cos(</v 2 - y> ri ) 

COS (ip r2 - (f Sa ). 



r 2 ,n 
-k 1 



(28) 



For such a system, it is possible to show ( Suppes et al. , 20121 Appendix) that a 
response x is selected when 



E 

s fl ,ri 
E 

s a ,r 2 
E 

r±,r 2 
E 

E 

T 2 ,S a 

E 

r 2 ,ri 



and 



k 1 

k 1 

s a ,r 2 

k 1 

r lt r 2 



ax, 

—ax, 

-a, 

ax, 

—ax, 

-a, 



a\/ 1 — x 2 



-aVT 
0, 



i 
I 

r 2 ,s a 
I 

r 2 ,ri 



—ay 1 - : 
a\/ 1 — x 2 
0, 



(29) 
(30) 
(31) 
(32) 
(33) 
(34) 

(35) 
(36) 
(37) 
(38) 
(39) 
(40) 



where a is a coupling strength parameter that determines how fast the solution 
converges to the phase differences given in (fl"3)) . We have similar equations for 
Ss(t), namely 

ifs- = uj - kf_ ri sin (<^ s _ - ip ri ) 
~kZ,r 2 sin(^ S7r - tpr 2 ) 
-Ks,n cos 0?** ^ <Pri) 
-k*. COS (tp s - -<Pr 2 ), 



(41) 
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ip ri = w - kf s _sin (ip ri - (p Sir ) 



rx 



rx,ri 0111 Vfrx ^r 2 J 
-kl uS _COs{tp ri - Lp s _) 

-Av 1)f . a cos(p ri -<Pr 3 ), (42) 

ip r2 = u Q - k® tTl sin (ip r2 - ip ri ) 
-k^ 2!S _sin{(pr 2 - f Sw ) 
-*v 2 ,n cos((^ r2 - <p ri ) 

-k 1 cos (ip r2 - ip S7r ), (43) 



and the equivalent equations for ([29l) - (|40|) . 

Equations (J^HJ) (SOJ) determine a given phase relation for a response x in a 
continuous interval, but they do not necessarily model the discrete case of a se- 
lecting X over Y. To do so, let us recall that because of the stochasticity of the 
initial conditions, if a particular value x is conditioned through couplings 
(|40|) , such response is selected according to a smearing distribution K{x\k Sa ). 
In fact, if / (y) is a simple reinforcement given by a Dirac delta function cen- 
tered in z, i.e. f(y) = S(y — z), then the response distribution r{x) becomes 



asymptotically ( Supped . 1959h 



r(x) = J k Sa (x\y)S(y - z)dy = k s (x\z), 

i.e. r(x) coincides with the smearing distribution at the point of reinforcement 
z. Thus, by using a reinforcement density given by a Dirac delta function 
centered at z, we are able to obtain the shape of the smearing distribution from 
the oscillator dynamics. Figure [1] shows the density histogram for a MATLAB 
simulation of the three oscillator model. Not surprisingly, the model's smearing 
distribution k s {x\z) fits well a Gaussian 

1 («-*) 2 

k s (x\z) = 



2tt 

where a = 0.05. We should mention that the estimation of the smearing distri- 
bution from a reinforcement schedule given by the Dirac delta function is not 
realistic for real psychological experiments. However, because we are running 
simulations without taking into consideration other behavioral and environmen- 
tal aspects, we are able to extract k s (x\y) from it. 

With the smearing distribution, we may now describe how a stochastic de- 
cision making processes may happen in the brain. Let X and Y be two of the 
possible responses for a certain behavioral experiment. In our model, response 
X would be associated with an oscillator, r% (t), and Y with the other, (t). If 
h = 2A 2 and J 2 = 0, then the response would clearly be I\. Since — 1 < x < 1, 
x given by 1161 we may state that x > corresponds to a response X and x < 
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Response 



Figure 1: Histogram density of 300 oscillator models reinforced using the dynamics given in 
equations (121 C — (I23P for y = 0.5. The histogram was obtained using 6300 points corresponding 
to 300 oscillators and 21 reinforcement trials per set of oscillators (trials 40 to 60). All 
parameters used were the same as those described in lSuppes et al. | J2012D . and K' was set to 
4, 500. The fitted Gaussian has fi = 0.51 and a = 0.05, with p < 10 -3 . 

to response Y. If we do this, when the reinforcement is for a value different 
from 1 or — 1, there is a chance, given by the smearing distribution, for each 
response. In other words, say a stimulus s a is conditioned to — 1 < z < 1. Then, 

P(X) = [ k Sa {x\z)dx 
Jo 

1 1 J z \ 
~ - + -erf —=- 

2 2 \V2a) 

and 

P(Y) = J k Sa (x\z)dx 
2 2 \V2aJ 

for \z\ , <J <C 1, where erf (x) is the error function. So, we can use the expression 

z = V^o-erF 1 (2P(R 1 ) - 1) 

to estimate the reinforcement value z leading to the desired probabilities of 
responses. Let us summarize the oscillator implementation of the Won/Lost 
version. We started with two stimulus, either Won or Lost, corresponding to 
the oscillators s a (t) and Sa{t). Once s a (t) or Sa(t) are activated, the dynamics 
of the system leads to a selection of either X or Y, corresponding to the decisions 
of betting or not betting on the next step of the game. 

We now model the "Disjunctive Version." Because the response Win or Lost 
is not known, we assume that both oscillators s a (t) and s-a (t) are activated 
simultaneously. Furthermore, since the subject also knows that we either Win 
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or we Lose, oscillators s a (t) and Sa (t) are incompatible, in the sense that they 
must represent off-phase oscillations in an SR-oscillator model. Therefore, the 
response system is not composed not only of three oscillators, but four: the two 
stimulus and the two response oscillators. The equivalent dynamical expressions 
for (f!?6")) - (|4!5)) when both oscillators are activated are the following. 

^sa = w o - kf a ri sin (tp Sa - ip ri ) 

-ks a ,r 2 sin (^a - <Pr 2 ) 
- k L,ri cos (^o ~ <Pn) 

-H a ,r 2 COs(ip Sa -<f r2 ), (44) 

(p s - = ujq- kf_ ri sin (ifs- - ip ri ) 

~kf_ T2 Sin ((Par-p ra ) 
-fcf_ Sa COS (<p s - - (f s J 
-fcf_ ri COS (lf s - ~ <p ri ) 

-fc^r 2 c °s(v s --^ 2 ), (45) 

ip ri = w - kf uSa sin (<p ri - ip Sa ) 
-kf u& _sax{ip ri - <p s -) 

-^n, s „ cos (ip ri - tp Sa ) 
-k I ri S _cos(L Pri - <^ s _) 

- fc n,r 2 C0S (^i ~<Pr 2 ), (46) 

ifir 2 = u a - fc^ ri sin (<p r2 - ip ri ) 

~ k r 2 ,s a Sin(<^r 2 - <ps a ) 

- kf 2 ^ sin ((f r2 - (f s -) 
-*4,n cos(^ r2 - ip ri ) 

- k r 2 ,s a COs((^ r2 - <p Sa ) 

-ki^cos^-ip^), (47) 

with the couplings between s a and r\ and ti as well as % and ri and r2 being 
the same, but 

k L^ = a '- («) 
As mentioned above, the parameter a' in (|48|) represents somehow a measure 
of the the degree of incompatibility of stimulus oscillators s a (t) and s„ (t) . 
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Now, as above, let I[ be the mean intensity at r\ (t) given by 

I'l = ((s Q (t) + S a(<) + ri(t)) 2 ), 

where this time we have both s a (t) and s-a (t) contributing to the intensity. 
Defining 8Lp ri a = ip ri — ip Sa and 5tp rij a = ip ri — Pa^, it is straightforward to 
show that 

I'l = A 2 ^+ COS (8tp ri :a ) + COS (5<f ru a) + COS (5ip ru a~ S(p ru a)^ ■ 

Similarly, for ri (t) we have 

4 = A 2 + cos (<fy?r 2 ,a) + COS (S(f r2i a) + COS {Sip r2 . a - 5<Pr 2 ,a) 



We notice that, because we now have more oscillators, extra interference terms 
shows up in the intensity. We now use the new intensities, /{ and I' 2 to compute 
the response, using 

b' = 



V - T' 



I'l + r* 

At this point, we should remark on the appearance of quantum-like effects. 
The model we are using, of inter fering o scillat ors, carry some similarities with 



the famous two-slit experiment (|Dirad . 119471 ). In the two-slit experiment, a 



quantum particle can go through two slits, and then hit a phosphorous screen, 
thus being detected. Since a particle is a localized object, going through one slit 
should mean no interaction with the other slit, and therefore we should expect 
to make no difference whether we keep the other slit closed or open. However, 
if we open both slits, we observe different probability patterns than if we close 
one of the slits. In other words, if we close one of the slits, and therefore know 
through which slit the particle went through, we get a probability of observing 
the particle in a certain region of the screen. If, on the other hand, we do not 
close one of the slits, we do not know through which slit it went, and when we 
do not know, the probability of detecting the particle changes. Formally, let L 
be the proposition "went through the slit on the left because the slit on the right 
was closed," R ="went through the slit on the right because the slit on the left 
was closed," and A ="was observed at a certain region on the screen." Then, it 
is possible to choose a region on the screen such that 

P{A\R) > P(^A\R) , 

and 

P(A\L)>P( r A\L), 

i.e., such that there are more particles being detected on it than outside of it. 
However, because of quantum interference, if we do not close the slits, and do 
not know where they went, it is again possible to have an A such that the above 
holds and 

P{A) <P (-.A) . 
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This is the two slit formal equivalent of the violation of STP. Now, let us go 
back to the oscillator model. For the the Win/Lost Version, a single stimulus 
oscillator is active at a time. Such oscillator produces a dynamics that leads to 
a certain balance between the responses X and Y. This corresponds to the only 
one slit open, when there is no interference from the other slit. However, for the 
Disjunctive Version, both oscillators are active, and their activity interferes with 
each other, in a way similar to the interference from the two slit experiment. The 
strength of interference, i.e. the degree of coherence between the two oscillators, 
is determine by the inhibitory coupling parameter k J s B _. Therefore, when both 
oscillators are active, we obtain a different probability distribution violating the 
standard axioms of probability, and consequently STP. 

To end this section, let us look at a simulation of the violation of the STP 
using neural oscillators. We simulated in Matlab R2012a the oscillators' dynam- 
ics determined by equations (|2l))) - P3")) and using the Dormand-Prince 
method. For simplicity, and without loss of generality, all oscillators were as- 
sumed to have the same frequency of 11 Hz, as this would allow the use of the 
couplings gi ven by (1551) - (PTO1) and (1151) without having to invoke learning equa- 
tions, as in ( Suppes et all 12012 ). For our system, the parameters utilized were 
the following: At r = 0.2 s, a v = y/njl, kf r = -kf = kf = k E 



<p V / ' S a ,ri S a ,V2 Tl,S a 7"2, s a 

—0.011 Hz, k ri T2 = k r2 ri = —11 Hz, k s _ ri — — k s _ r2 — — k ri s _ — k r2 S _ = 1, 

h E — —h E — h E — h E — — fl ("117 H7 h 1 — —h 1 — —h 1 — 

H 2 ,s w = 1) H t ,r 2 = K 2 ,r! = 0, and kl a S _ = fcf_ Sa = -0.011, where At r is 
the time of response after the onset of stimulus, and a v is the variance of the 
initial conditions for the phase oscillators. We ran 1,000 trials with the above 
parameter values, for the Win/Lost Version we obtained X as response 63% 
when the Lost oscillator was used and 58% with the Win oscillator. However, 
for the Disjunctive Version, when the two oscillators were active, we obtained 
X as a response 36% of the runs, clearly showing a violation of Savage's STP 
and an interference effect. 



5. Conclusions 

In this paper we presented a neural oscillator model based on reasonable 
assumptions about the behavior of coupled neurons. We then showed that 
the predictions of such model presents interference effects between incompat- 
ible stimulus oscillators. Such interference leads to computation of responses 
which are inconsistent with standard probability laws, but are cons i stent with 
quantum- like effects in decision making processes Busemever et al. ( 20091 ). Of 



course, the fact that we get interference from si and S2 does not preclude us from 
necessarily making rational decisions. We could, for example, modify the above 
model, and include an extra oscillator, such that after reinforcement we could 
eliminate interference, therefore satisfying the STP. But our model suggests that 
a plausible mechanism for quantum effects i n the brain ma y not need to rely on 
actual quantum processes, as advocated by IPenrose] (1989), but instead can be 
a consequence of neuronal "interference." 
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One interesting feature of the oscillator model is that the amount of inter- 
ference between the oscillators is determined by the couplings between si and 
s-2- The stronger the incompatibility between the two stimulus, the more we 
should see interference between the responses. Furthermore, because interfer- 
ence is a consequence of inhibitory couplings, we could in principle design exper- 
iments where pharmacological interventions could suppress inhibitory synapses, 
destroying such interference effects. 

It is worthwhile to compare our model to the one proposed by iKhrennikov 
(2011). Whereas in Khrennikov's model the interference of electromagnetic 
fields account for quantum-like effects in the brain, in our model we use Ku- 
ramoto oscillators to describe the evolution of coupled sets of dynamical os- 
cillators close to Andronov-Hopf bifurcations. Because such oscillators obey a 
wave equation, they are subject to interference. So, in some sense, we may 
say that the oscillator and the field models are equivalent representations of 
quantum-like interference in the brain. However, though we can think of the 
weak couplings between oscillators originating from the synapses between neu- 
rons, it is also conceivable that the weak oscillating electric fields generated 
by sets of oscillators provide sufficient additional coupling to account for their 
synchronization and interference. Thus, we believe that there are not equiva- 
lencies between the models, but also that our oscillator model might provide an 
underlying neurophysiological justification for Khrennikov's quantum-like field 
processing. However, we point out that because Khrennikov's model relies on 
a quantum representation of states in terms of Hilbert spaces, his model brings 
a more powerful machinery. It would be interesting to prove some represen- 
ta tion theorem bet ween neural oscillators and the field model, but, as shown 
m Ide Barrod (|2012h . this seems implausible given that oscillators may generate 
processes that are not representable in terms of observables in a Hilbert space. 
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